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FOREWORD 


This  report  is  one  of  a  series  of  related  papers  covering  various 
aspects  of  a  broad  program  to  investigate  the  flow-  field  variables 
associated  with  hypersonic -velocity  projectiles  in  free  flight  under 
controlled  environmental  conditions,  Ine  experimental  research 
is  being  conducted  in  the  Flight  Physics  Range  of  GM  Defense 
Research  Laboratories,  General  Motors  Corporation,  and  is  sup¬ 
ported  by  the  Advanced  Research  Projects  Agency  under  Contract 
No.  DA-dl-021-AMC-11358iZ).  It  is  intended  that  this  series  of 
reports,  when  completed,  will  provide  a  background  of  knowledge 
of  the  phenomena  involved  in  the  basic  study  and  thus  aid  in  a 
better  understanding  of  the  data  obtained  in  the  investigation. 
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ABSTRACT 


An  exact  solution  to  the  heat  conduction  equation  is  developed  for  a 
solid  spherical  region  subjected  to  an  arbitrary  time -independent  surface 
heat  flux.  It  is  assumed  that  the  thermal  properties  are  independent  of 
temperature  and  that  no  sources  or  sinks  exist  in  the  region.  It  is  shown 
that  the  solution  is  valid  for  heated,  cooled,  or  thermally  insulated  regions. 
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NOMENCLATURE 


SUBSCRIPTS 

o 

i 

n,  m,  J 


radius  of  the  sphere 

expansion  coefficients 

specific  heat 
given  constant 

angular  dependence  of  heat  flux 
thermal  conductivity 
constant 

associated  Legendre  function 

heat  flux  at  surface  of  sphere 

heat  generation  per  unit  volume  per  unit  time 

nondimensional  radius 

nondimens  tonal  time 

temperature 

cos  8 

space  variable 

thermal  diffusivity 

constant,  defined  by  Equation  (7) 

density 

spherical  coordinates 
eigenvalue 

time 

initial  conditions 
space  coordinate 
expansion  indices 
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INTRODUCTION 


Transient  heat  conduction  problems  in  polar  spheric  oordinate  systems 
are  of  considerable  importance.  One  important  class  of  boundary  conditions  for 
the  spherical  region  is  a  prescribed  heat  flux  at  the  surface.  Since  the  heat 
flux  during  the  initial  phases  of  reentry  is  approximately  independent  of  time, 
the  solution  to  this  class  of  problems  has  application  to  the  transient  heating 
regime  of  reentry  vehicles.  In  particular,  the  solution  of  the  arbitrary  surface 
flux  problem  can  be  used  to  calculate  the  time  for  incipient  phase  change  of 
spherical  models  utilized  to  "^ulate  reentry  in  free -flight  range  experiments. 
The  temperature  distribution  in  a  spherical  region  subjected  to  a  constant  heat 
flux  has  been  reported  in  References  1  and  2.  However,  the  problem  of  a 
spherical  region  subjected  to  a  variable  surface  flux  has  not  been  investigated. 

Considerable  research  has  been  conducted  concerning  transient  tempera¬ 
ture  distributions  in  bodies  undergoing  phase  changes.  References  2  through  9, 
including  the  references  therein,  provide  a  relatively  oomprehetisive  biblio¬ 
graphy  of  this  subject.  However,  it  is  to  be  noted  that  the  above  authors  assumed 
that  the  transient  heat  conduction  could  be  represented  by  a  one-dimensional 
model.  When  the  heat  flux  is  a  function  of  position,  various  sectors  of  the  sur¬ 
face  will  reach  the  phase  transformation  state  at  different  times.  Hence,  the 
phase  transformation  is  further  complicated  by  this  effect,  and  it  is  important 
to  be  able  to  evaluate  the  surface  temperature  distribution  and  time  when  this 
temperature  is  first  reached  at  some  point  on  the  body. 

The  purpose  of  the  present  paper  is  to  determine  the  exact  transient  tempera¬ 
ture  distribution  for  a  spherical  region  subjected  to  a  time-independent  surface 
heat  flux  that  may  be  a  function  of  position.  The  solution  's  obtained  for  a 
material  with  constant  thermal  properties  and  a  region  without  sources  or  sinks. 
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FORMULATION  OF  THF  PROBLEM 


The  partial  differential  equation  for  feat  conduction  in  an  isotropic,  solid, 
including  the  effects  of  heat  generation  anrf  variable  thermal  properties,  can 
be  written  as 

PC  =  v  (KvTltQI^  ,  r)  (1) 

The  coordinate  system  for  this  problem  is  assumed  to  be  the  region  bounded  by 
the  sphere  of  radius  a  as  shown  in  Figure  1. 


x 

Figure  1  Coordinate  System 

In  order  to  simplify  the  solution  of  Equation  (1)  it  in  assumed  that:  (a)  the 
heat  generation  per  unit  volume  per  unit  time  is  aero;  (b)  the  thermal  properties 
are  constant;  and  (c)  the  boundary  conditions  are  independent  of  time. 

It  is  assumed  that  ihe  initial  temperature  distribution  is  known  and  that 
an  arbitrary  heat  flux  at  the  surface  of  the  sphere  is  specified.  Utt  zlttg  the 
three  assumptions  listed  above,  Equation  (1)  can  be  written  as 


(2) 


m  m  eirsKss  seskarck  l  a  a  ©  a  a  t©  is  5  c  s  ®  oixetAi  motors  co^r&RArtOR 


mm-mu 


The  solution  of  Equation  (2)  is  subjected  to  the  following  initial  and  boundary 

f'lin^HnnK 


t (c »  e,  0,  o)  =  tq(C,  e,  0) 


-(-*£) 


^'c-a 


=  Fg(e,  0) 


In  Equation  (3b),  F  is  known  constant  and  K  is  the  thermal  conductivity  of 
the  material  In  order  to  generalize  the  solution  of  Equation  (2)  subjected  to 
(3a)  and  (3b),  we  define  the  following  nondimensional  distance  and  time 

r  =  C/a  (4 


Combining  Equations  (2),  (4a)  and  (4b)  and  expanding  the  Laplaclan  in  terms  of 
spherical  coordinates  yields 


BT  B2T  2  ST  1  S  /  .  .  BT\ 


sin  6  50 


5?T 

re¬ 


converting  the  boundary  conditions  to  tee  nondlmens ionoiized  form  gives 


T  (r,  0,  0,  o)  -  Tq  (r,  0,  0) 


0  g{e,0) 


where 


The  solution  of  Equation  (5)  must  also  most  the  physical  const*  ctnt  that  the 
temperature  throughout  the  region  is  bounded  for  finite  time. 
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SOLUTION  OF  THE  PROBLEM 


It  san  be  shown  that  a  simple  variables-separable  solution  of  Equation  (5) 
cannot  meet  all  of  the  boundary  conditions.  This  implies  tha'.  if  a  solution  with 
separated  variables  is  valid  for  die  given  system,  the":  re  must  be  a  composite 
function. 

Assume  that  the  solution  of  Equation  (5)  can  be  written  in  the  following 
form 


T(r,  e,  0,  t)  =  f(t)  +  u(r  ,  8,  0)  +  v  (r,  8,  0,  t)  . 


(8) 


It  can  be  s^own  that  a  solution  to  Equation  (5)  that  meets  (6a)  and  (fib)  lb  unique. 
A  mathematical -uniquene&a  proof  for  this  system  is?  given  by  Churchill  The 
physical  uniqueness  of  the  solution  can  also  be  proven  by  invoking  the  Second 
Law  of  Thermodynamics.  Assume  that  Tj  and  Tg  are  unique  solutions  to 
Equation  (5)  subjected  to  (8a)  and  (fib).  If  T-  and  are  solutions,  then 
(TrT?)  mus  also  be  a  solution.  It  can  be  easily  shown  that  the  Second  Law 
of  Thermodynamics  is  violated  if  the  quantity  Tj  -  ?g)  i«  not  identically  zero. 

Combining  Equations  (5)  and  (8)  yields 

f  (t)  +  |^  =  v2u  +  V2 V  .  (9) 


Let 


I 


f’(t)  =  '2u 


(10) 


dv 

¥T  = 


V2  V  . 


(li) 
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Since  the  right-hand  side  of  Equation  (10)  is  independent  of  time,  the  left-hand 
side  must  also  be  independent  of  time.  Therefore,  i(t)  can  be  written  as 

f(t)  =  Mt  (12) 

where  M  is  a  constant.  The  boundary  and  initial  conditions  of  the  problem  are 
satisfied  by 


^  (r,  6,  0)  =  u(r,  0,  0)  + 

H.t  ‘  ?g<8,0)  '  5?l 

Equations  (1S^  and  (14)  are  satisfied  if 

If  I  •**«.») 

'r  =  l 


v(r,  8,  0,  o) 


r  =  l 


By 

Br 


'r-i 


we  arbitrarily  set 


By 

Br 


=  0 


and 


(13) 

(14) 


(15a) 


(15b) 


v  (r,  6,  0,  o)  =  Tq  -u(r,  8,  0)  .  (16) 

Therefore,  the  solution  to  the  initial  problem  is  reduced  to  the  solution  of 
the  following  two  differential  equations  and  their  corresponding  initial  and 
boundary  conditions: 


Br 

at'  “ 


y2  v 


(17a) 


subjected  to 


By 

Br 


r  =  l 


0 


(17b) 
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v(r,  6,  0,  o)  =  TQ  -u(r,  g,  0  } 
and 

v2  u  =  M 
subjected  to 

I  =  P  g(  6,  0)  .  (18b) 

ar  'r  =  1 

A  complete  solution  co  the  initial  problem  is  obtained  when  the  solutions  u  and 
v  are  evaluated  and  the  constant  M  in  Equation  (12)  is  determined. 


(17c) 

(18a) 


SOLUTION  OF  POISSON'S  EQUATION 

The  solution  of  Equation  (18a)  c>  .  he  obtained  by  making  the  following 
definition: 

u(r,  g,  0)  =  w(r,  9,  0)  +  «<(r)  (19) 

where 

v2  to  =  0  (20) 

a>.d 

+  jp  'Jr  *  M  .  (21) 

Equation  (21)  is  a  second-order  nonhomogeneous  differ  ential  equation  whose 
general  solution  can  be  written  as 

♦fr)  “  r  ^  +  ‘7~  +  ^2  *  (^2) 

Since  u  must  remain  finite  at  the  origin,  the  value  of  must  be  identically 
zero.  Arbitrarily  setting  the  value  of  $  equal  to  zero  at  the  surface  of  the 
sphere,  we  obtain 

*(r)  =  ^  (r2  -  1)  .  (23) 
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.The  solution  of  th  w  equation  can  be  obtained  by  separation  of  variables 
w(r,  0,0)  =  R(r)6'  )*0)  (U) 


Combining  Equations  (20)  and  (24)  and  dividing  by  u>  yields 


2  R"  2rR 
r  IT- +  ~R~ 


(l~x2)  +  (k*J 


9  (  (\  xV6\-  -  *' 
^ ( (1  x >35 )  5 


m 


(isb) 


where 

x  =  cos  0  (28) 

and 

m  =  a  separation  constant. 

Therefore,  the  solution  for  the  4>  function  can  be  written  as 

*»)  =  (27> 

Equation  (25)  can  be  rewritten  as 

r2R"  +  2rR  +  J32  R  =  0  (28) 


and 

4-  [<‘-x2)  M] +  [^ '  ] e  “ 0  (29) 

where 


/J2  is  a  separation  constant. 

The  solution  to  Equation  (29)  Is  given  by  the  associated  Legendre  function 


©(x)  =  P“  (x) 


(20) 


and 


$  =  n(n+l) 


n  =  0,  1,  2  — 


(31) 
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Equation  (38)  is  Cauchy's  linear  differential  equation  which  has  the  solution 

Rn  (r)  =  C?  rn  +  C4  r  "(n+1),  (32) 

For  the  solution  to  be  bounded  at  the  origin,  must  be  identically  Eero. 
Hence  the  solution  to  Equation  (19)  can  be  written  as 


u(r,  9,  0)  =  jp(r2-l)  +  £  £  yI1  Pj^OO  |AnmC0S  +Bnmsin  ^ 

n=o  ni=o  (  j 

Applying  boundary  condition  (15a)  to  Equation  (31)  yields 

oo  n  /  i 

0g(8,  0)  =  S  nPn^  iAnmC0S  +  BnmSlnm0^  ’  ^ 

n=l  m=o  (  ) 

Now  the  arbitrary  function  (0  g (9,  0))  can  be  expanded  in  the  following 
manner: 

0g(M)=£  £  P“(x)  {a  cosra0^bnmslnm0|.  (35) 

n=o  m=o  (  ; 

The  numerical  values  of  the  expansion  coefficients  a _ and  b _ can  be 

nm  nm 

obtained  by  invoking  the  orthogonality  relations  between  the  trigonometric  and 
associated  Legendre  functions. 


The  normalization  integrals  for  the  associated  Legendre  functions  are 
given  by 


2  (n+m)  ! 

(1 £  n+I}  (n-rn)! 


(36) 


and  the  values  of 

nm 


and  b 


nm 


are  given  by  Equations  (37)  Oirough  (36). 
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1  2V 

aoo  *hff  $  S*  fa  t)  d*  d  *  (37) 

-1  *0 

anm  =  %rrEm)"T)  1  /  Pn  (x)  /  0)  ccsm0  dp  dx  (38) 

-1  o 

1  2? 

bnm  =  ^THWr1  /  P“(x)  / 0g(^0)slnm0d0dx  .  (38) 

-1  o 

The  vaiues  of  the  unknown  coefficients  in  Equation  (34)  are  obtained  by  equating 
like  terms  in  Equations  (34)  and  (35).  Therefore,  it  is  seen  that 


M  = 

3a 

00 

(40) 

Anm 

an m 
n 

n  4* 

0 

(41) 

Bnm 

b 

nm 

n 

n  =£ 

0  • 

(42) 

Since  the  values  of  a„_  and  b _ are  determined  in  terms  of  the  given 

run  nm  9 

normal  boundary  condition  and  simple  quadratures,  the  solution  of  the  u  function 
is  determined.  It  is  also  to  be  noted  that  the  constant  M  (see  Equations  (8)  and 
(12) )  has  been  specified  in  terms  of  the  normal  bounds*  y  condition. 
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SOLUTION  -Oh'  THE  v  EQUATION 

The  purpose  of  this  section  is  to  outline  the  solution  of  Equation  (17a) 
subject  to  the  boundary  and  initial  conditions  as  given  by  Equations  (17b)  and 
(17c).  The  solution  is  obtained  by  assuming  that  the  function  is  separable  and 
by  proceeding  in  the  usual  manner.  This  attack  yields  the  following  results: 

v(r,  0,  0,  t)  =  R(r)  6(0)  #(0)  e  X  (43) 

where 


Equation  (48)  is  similar  to  Bessel's  equation,  and  it  can  be  put  into  the  correal 
form  by  the  following  change  of  variables. 

Let, 

R(r)  =  ^172  '  (47) 

Combining  Equations  (46)  and  (47)  yields 


10 
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The  solution  to  equation  (48)  is  given  by  the  Bessel  (unction  of  order  (njT/2). 


G  (r)  =  J 


n+1/2 


<X„r) 


(49) 


Combining  (47)  and  (49)  yields 


R(r) 


Jn+l/2  *Xnr* 

=  T7T- 


(50) 


Since  the  solution  must  meet  die  boundary  condition  (17b)  which  is  independent 
of  0,  0,  and  t,  the  derivative  of  R(r)  evaluated  at  r  =■- 1  must  be  zero. 


r'(i) 


d 

3r 


r 

1? 

*5* 

<M 

v-4 

-* 

\ . 

L  r*'*  J 

(51) 


It  can  be  shown  that  Equation  (51)  is  satisfied  for  the  eigenvalues, 
given  by 


nj 


t 


nJn+l/2*\ij*  =  Jn+3/2  ( *hj } 


(52) 


where 


n  •*  0,  1,  2  — 

J  =  1,  2  —  . 


The  first  ten  eigenvalues  X  are  given  for  n  equals  one  to  fifteen  in  the 
appendix  to  this  report.  Therefore,  the  solution  to  Equation  (17a)  that  meets 
the  normal  boundary  condition  (17b)  is  given  by 


cd  oo  n  J  .,,(1  .r) 


CnmJ  C0S  +  DnmJ  Bhl  m  0 


’4 


(53) 


li 
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The  numerical  value  of  the  coefficients  C^,  and  can  be  evaluated 


from  the  initial  conditions  as  given  by  Equation  (17c). 

(T„-u) .  g  £  £ 

n=o  j=l  m*o  r  '  ' 


cos  m0  +  sinm0 1 


Since  sin  m  0  ,  cos  m  0 ,  and  p“(x)  form  orthogonal  sets  of  functions  and 
the  Bessel  functions  are  orthogonal  with  respect  to  the  weight  function  r ,  the 
values  of  CRmj  and  can  be  evaluated  in  terms  of  quadratures  It  can 

be  shown  that  the  value  of  the  normalization  integral  for  the  fractional-order 
Bessel  functions  subjected  to  the  eigenvalues  of  Equation  (52)  are  given  by 


/r[ WV!]2<lr=  “ 


n(n+1)  Vl/a<V 


The  values  of  and  Dnmj  are  determined  to  be 


fC-l) 


X*  (2n+l)  (n-m) ! 

«|" 

J  VOl 

J  f 1  Ji  ff 

/rV2  Vl/2-'^r)  /  <x> /<V">  joS  j  aad,[dr  • 


The  constant  n  must  be  replaced  by  2v  when  is  evaluated. 
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RESULT 

The  solution  fs  fi  ’ally  obtained  by  combining  the  various  parts. 


T-Mt  +  $  (r2-l)  +  f;  irnp”(s)  (a 
n~o  m-o  ( 


nan  co®  m3  +  Bnm8iKm0 


co  n  oo  J  wo(X^,r)  (  i  -f 

+  EEE  "+  VI  Pn  (x)  I  CnmJ cos  m0*Dnmj  alam0l  0  *  ** 

n-o  m=o  j=i  r  '  1 


-x2_,t 


(5?) 


The  values  of  the  expansion  coefficients  M.  A  ,  B  ,  C  .  D 

nra’  nm  ’  nmj  •  nmj 
and  the  eigenvalues  are  given  by  Equations  (40),  (41),  (42),  (58),  and  (52), 


The  results  of  Equation  (57)  can  be  simplified  if  the  initial  temperature 
distribution  is  assumed  constant. 


In  this  case,  the  left-hand  side  of  Equation  (57)  can  be  replaced  by  the  tempera¬ 
ture  excess  (T-Tq),  and  the  function  (Tq-u)  in  the  integrands  of  the  expansion 
coefficients  C^.  and  Dnmj  can  be  replaced  by  (-u). 

If  the  heat  flux  and  temperature  distribution  are  axially  symmetric  functions, 
the  solution  and  boundary  conditions  are  not  dependent  on  the  angle  0 .  This 
simplification  1  reflected  by  setting  the  value  of  m  identically  equal  to  zero 
in  Equation  (5».  J  in  addition  to  the  above  the  initial  temperature  of  the  region 
is  constant,  further  simplifications  result  in  the  solution. 


@M  SiriNI£  ftKSSAftCH  LA#ORATOnig$  $  SESKfJAL  a6TO*S  C  ©  tf  #0 *  A 1 1  ©  H 


TR85-01M 


DISCUSSION  OF  THE  RESULTS 


Some  comments  may  be  made  concerning  the  generality  of  the  above  solution. 

The  development  in  the  proceeding  sections  has  shown  that  the  temperature 
distribution  in  a  spherical  region  subjected  to  a  known  heat  flux  at  the  boundary 
can  be  represented  in  terms  of  known  functions  and  quadratures.  It  is  also 
important  to  note  that  the  eigenvalues,  X^j  ,  of  Equation  (52)  are  not  dependent 
on  the  form  of  the  given  normal  heat  flux  distribution.  Therefore,  once  the  set 
of  eigenvalues  are  determined,  they  are  valid  for  all  input  functions  g(  0,  0). 

The  generality  of  the  above  solution  is  manifested  by  the  expansion  coefficients 

a  and  b„m  .  The  sign  and  numerical  value  of  a  indicates  the  character 
nm  run  oo 

of  the  problem  that  is  being  studied.  If  the  value  of  aQQ  is  greater  than  zero, 
there  is  a  net  flux  of  energy  into  the  region  and  one  would  expect  the  internal 
energy  of  the  region  to  rise  with  time.  This  is  confirmed  by  the  first  term  in 
the  solution.  On  the  other  hand,  the  numerical  value  of  a  is  negative  if  there 

OO 

is  a  net  flux  of  energy  leaving  the  region.  Physically,  the  negative  sign  implies 

that  the  internal  energy  of  the  region  decreases  with  time.  For  a  system  in 

which  the  integral  of  the  flux  over  the  boundary  is  zero,  one  expects  that  the 

net  internal  energy  of  the  region  wouia  remain  constant.  This  phenomenon  is 

also  predicted  by  the  solution,  since  under  the  above  assumption  the  value  of 

would  be  zero.  Finally,  If  the  heat  flux  is  identically  zero  over  the  whole 

surface,  all  of  the  a  's  and  b  's  are  zero,  and  the  solution  reduces  to 
am  nm 

the  temperature  history  of  an  insulated  spherical  region  with  an  initial  tempera¬ 
ture  distribution^ 

The  above  discussion  implies  that  the  temperature  variation  given  by  Equation (57) 
is  valid  for  a  wide  range  of  normal  boundary  conditions  that  include  net  heating, 
cooling,  or  an  insulated  boundary. 
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APPENDIX 

FIRST TFN  EIGENVALUES  FOR 
n  EQUALS  ONE  TO  FIFTEEN 


10  FIRST  JO  EIGENVALUES  FOR  »  Oi 1 » 2 • . * .  .  1  5 
N*  0  T  '+4493432*5+0 1  +77252684+01  4  1 0904  1 99+02+ 14066 1  1  3+02+ 1  7220751+02“ 

+20371 386+02+23519482+02+26666 2 1 0+02+2981 1 572+02+3295634 7+02 
ft%\  +208 l 5673+0 1 +59403 >4  O l  +  92058377+0 1  + 1 2404468+02+ 15579242+02 

+ 1 8742647+02+2 1 090700+02+280 S2822J  02  +  28203759+02+31 352090+02 

+33421  786+0* +728' >9351+  >1  +1061  386!  +02+138*61  09+02 +T 7042904  +02 
+  2022 1  359+02+2.3790484+02+265526!  2+('2+2971  0290+02  +  .VR3AH32+0? 

#•3  +45  1  40868+0  1  +85837801  +C  1+1  I  97274  t  +02+  1  524431  6+02-rli  8468 1  48+02 

<•21666604+07+24850079+02+28023875+02+31 191211 +02+34353380+02 
t(* «  "  +564669 18+01 +98404477+0 1 + 1 3295565+02+1 6609347+02+ 1 98624 14+02 ' 

+23082797+02+26283269+02+29470638+02+3264889 1 +02+35J.2054 1 +02 

__  +67565428+0  1  +  1 1 07020R+02  +  1 4590564+02+T  7947J 73+02+2 1 231 072+0? 
+24474R26+ 07+2 76977 19+02+30895998+02+ 74 086597+0?+372686?8+02_ 

“  *»6  +7851 1229+0  1  +  1 2279339+02+ 1586.7226+02+1926271 3+02+22376/ 60+02 

+258460R2+02+290 ’4347+02+32302513+02+35506331  02+38699590+02 
it*  7  +8974  8507+0 1  + 1 3472027+02+ 1711751 6+02+20559426+02+23906448+02 

+27199. 55+OP+30487507+02+33692! 73+02+36909921 +02+401 1 507! +02 
**  8  '  "  + 1 OO 10363+02+ 1 4651260+02+ 1 8356320+02+21 840002+02+25218701 +02 

+28536456+02+31815102+02+35066761+02+3829891 7+02+41 51 6645+"7 
#“9  +11 079406+02+1 58 1 9226+02+ i 9581887+02+23 1 06579+02+26516699+02 

+29R59493+02+33 1 58755+02+364277 1 6+02+3967462 1 +02+4290*926+02 
“  <*10  +"l  2  1 431 88+02+T6977545+02+26795968+02+24360766+02+2 7801 882+02 

+  31 169808+02+3448981 8+02+37776220+02+4  I  038202+02+44261 551+02 
*11  ”  +1 3202624+02+1 S127575+02+?iO99954+d2+25604 055+02 +290 758! 1+02 

+32468636+0?+35809372v02+391 1 3474+02+42390605+02+45647219+02 
Ate!  2  + 14258328+02+ 19270297+02+23 195006+02 +2683751 5+02+3C339.26+02 

I  33757027+02+371  18471  +  02+4O440269+02+43732665+02+47002744  +  02 
~ tfs  I  3  ""  _  +153 10864 +02+20406585+02 +2*^8204 0+02+280621  45+^02+3 1 593688+02 

+3503587t,+ 02+384 1 7952+02+41757424+02+45065265+02+48349022+02 
Ate  1  4  +  *  6360680+02+2 1537115+02  +2556 1  d  M+02+29..  78732+02+32839777+02 

+363C 5993+ 02+ 39708552+02+4 3065683+ 02 +4 6388927+02+496861 05+02 
it*  15  +T?4 08 044 +02+226624 9 w+02 +26735 168+02 +30488006+ 02+34077861 +02 

+37568007+02+40990923+02+44365655+02+4774  *39+02+51 G 1 31 84+02 


X  given  for  n  -  0,  1,  2,  ....  15  q  =  1,  2,  3,  . . .  10 

it,  q 


where: 


X  are  the  set  of  eigenvalues  >  0  to  the  following  equation. 
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